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Grant  No:  AFOSR-90-020 1 

Title:  Dynamic  modeling  of  joint  lifednies,  of  marginalization  and  of  multi-unit  minimal  rc  - 

pairs. 

Due  date:  30  June  1993 

The  proposal  for  Grant  AFOSR-90-0201  consisted  of  the  following  parts: 
x-  (1)  An  introduction  and  a  study  of  multivariate  mean  residual  life  functions. 

i  (2)  A  study  of  some  dynamic  notions  of  multivariate  aging. 

(3)  Modeling  dynamic  conjunction,  separation,  and  marginalization  of  reliability  systems. 

(4)  A  further  study  of  the  multivariate  notions  of  minimal  repairs. 

i  We  have  achieved  most  of  the  objectives  of  Utc  proposal  during  the  last  three  years. 

L 

(1)  Consider  a  multi-unit  system  consisting  of  a  number  of  components.  In  order  to  predict  and 
study  the  stochastic  behaviour  of  the  system  one  usually  need  to  model  the  joint  distribution  of 
the  lifetimes  of  the  components.  There  are  several  ways  in  which  tliis  can  be  done.  The  clas¬ 
sical  approach  is  to  simply  postulate  a  “leasonable”  multivariate  joint  distribution  and  hope  that 
the  postulated  distribution  indeed  corresponds  closely  to  the  real  life  situation  which  it  tries  to 
describe.  Another  common  approach  is  to  postulate  the  Laplace  transform  of  the  “reasonable” 
life  distribution,  Our  approach  in  the  last  few  years  has  been  to  postulate  the  underlying  life 
distribution  through  the  intuitive  notion  of  the  multivariate  conditional  hazard  rate  functions. 
We  have  shown  that  the  multivariate  conditional  hazard  rate  functions  uniquely  determine  the 
joint  distribudon  of  the  underlying  lifetimes  and  this  justifies  the  use  of  die  multivariate  condi¬ 
tional  hazard  rate  functions  for  the  purpose  of  modeling.  We  recently  have  suspected  that  a 
similar  modeling  procedure  can  be  introduced  by  the  use  of  the  multivariate  conditional  mean 
residual  life  functions  rather  than  the  multivariate  conditional  hazard  rate  functions.  It  turns 
out  that  indeed  this  is  possible  in  many  situations,  tiiough  with  some  exceptions.  The  results 
of  our  investigations  in  diis  area  arc  described  in  [1],  A  critical  comparison  of  these  two  dis¬ 
tinct,  though  similar,  approaches  can  be  found  in  [2].  We  have  also  continued  our  study  of 
the  multivariate  conditional  hazard  rate  functions.  The  use  of  the  multivariate  conditional  haz¬ 
ard  rate  functions  for  the  purpose  of  simulation  is  described  in  [3].  A  discrete  version  of  these 
functions  is  introduced  and  studied  in  [4J.  The  latter  paper  is  based  on  a  pan  of  tlie  Ph.  D. 
dissertation  of  the  doctoral  student,  Jose  Valdcz-Torres,  supervised  by  one  of  the  Pi’s  (Moshc 
Shaked).  In  [5]  we  introduce  and  study  a  model  of  failures  of  units  with  discrete  lifetimes, 
We  suppose  that  a  unit  has  a  sequence  of  ta.sks  to  perform  and  that  it.s  lifetime  is  measured  by 
the  number  of  tasks  performed  before  its  final,  fatal  failure.  Upon  a  failure  the  unit  may  be  re¬ 
paired  (with  some  probability)  and  then  it  is  given  another  chance  to  perform  the  current  task. 
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The  unit  dies  when  (with  some  probability)  a  repair  cannot  be  completed.  We  derive  some 
stochastic  comparisons  of  pairs  of  such  models.  Tlie  stochastic  comparisons  are  then  applied 
for  obtaining  results  regarding  the  inheritance  of  several  aging  properties  by  the  repaired  unit. 
Various  examples  illustrate  the  applicability  of  the  model.  Some  variants  of  the  model  of  this 
paper  can  be  viewed  as  discrete  analogs  of  the  notion  of  imperfect  repair.  The  discrete  multi¬ 
variate  hazard  rate  functions  (which  are  closely  related  to  the  multivariate  mean  residual  life 
functions)  are  the  basic  technical  tool  that  is  used  in  this  paper,  Anotlier  work  in  this  area  is 
done  in  [6].  Consider  a  reliability  system  consisting  of  n  components.  Suppose  that  the  fail¬ 
ures  and  the  repair  completions  of  the  components  can  occur  only  at  positive  integer-valued 
times.  At  any  time  k  each  component  can  be  in  one  of  two  states:  up  (i.e.,  working)  or  down 
(i.e.,  failed  and  in  repair).  The  system  state  is  also  either  up  or  down  and  it  depends  on  the 
states  of  the  components  through  a  coherent  structiue  function.  In  this  paper  we  formulate 
mathematically  the  above  model  and  we  derive  some  of  its  properties.  In  particular,  we  identi¬ 
fy  conditions  under  which  the  first  failure  times  of  two  such  systems  can  be  stochastically  or¬ 
dered.  A  variety  of  special  cases  is  used  in  order  to  illustrate  the  applications  of  the  derived 
properties  of  the  model.  Some  instances  in  which  the  times  of  first  failure  have  the  NBU 
(=new  better  than  used)  property  are  pointed  out.  This  work  is  based  on  a  part  of  the  Ph.  D.. 
dissertation  of  the  doctoral  student,  Haolong  Zhu,  supervised  by  one  of  the  Pi’s  (Moshe 
Shaked). 

(2)  Most  of  the  univariate  notions  of  aging  such  as  logconcavity,  IFR,  NBU  etc.  can  be  character¬ 
ized  by  some  kind  of  a  univariate  stochastic  monotonicity  of  tlie  mean  residual  lives  of  the  cor¬ 
responding  underlying  components.  We  have  recently  have  developed  several  dynamic  no¬ 
tions  of  multivariate  stochastic  orders.  Having  such  orders  it  is  natural  to  try  to  apply  them  for 
tlie  purpose  of  developing  new  useful  and  intuitive  notions  of  multivariate  aging.  This  has 
been  done  in  [1 1].  A  discrete  version  of  these  multivaiiate  stochastic  orders  has  been  devel¬ 
oped  in  [8].  The  latter  paper  is  based  on  another  part  of  the  Ph.  D.  dissertation  of  Jose 
Valdez-Toires.  Having  lliese  discrete  multivariate  stochastic  notions  it  is  now  pos.sible  to  in¬ 
troduce  and  study  discrete  analogs  of  the  multivariate  aging  notions  which  are  studied  in  [5J. 
This  was  done  by  a  doctoral  student,  Jose  Rocha,  at  the  University  of  Arizona,  under  the  su- 
pei-vision  of  one  of  the  PI’^  (Moshe  Shaked).  Another  related  study  about  characterization  of 
aging  notions  through  block  replacement  policies  is  described  in  [9].  The  latter  paper  is  based 
on  another  part  of  the  Ph.  D.  dissertation  of  Haolong  Zhu.  In  [5]  these  notions  play  an  impor¬ 
tant  role  in  the  quest  for  understanding  the  aging  propcities  of  the  distributions  of  the  actual 
lifetimes  of  the  repaired  components.  In  [lOJ  we  obtain  stochastic  comparisons  of  some  quite 
general  wear  processes.  Let  {N^(0,^0)  and  {Ng(r),<^}  be  two  counting  processes  with 
respective  independent  intcrarrivals  Aj  and  ,  and  let  ( and  {Y^}  be  two  sequences  of  in¬ 
dependent  non-negative  random  variables.  Consider  the  two  compound  processes  which  are 
based  on  tliese  variables.  Then  it  is  well  known  iltat,  under  .some  assumptions,  if  the  Aj-  are  in- 
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dependent  of  iheX|,  and  if  the  B/  are  independent  of  the  Yp  then  the  resulting  processes  can 
be  compared  stochastically.  In  this  paper  we  provide  conditions  under  which  this  stochastic 
order  relation  holds  even  under  weaker  conditions.  In  fact,  we  derive  the  stochastic  order  re¬ 
lationships  among  processes  which  arc  much  more  general  titan  the  cumulative  ones  described 
above.  For  example,  it  is  not  necessarily  assumed  that  the  A,-  are  independent  of  theX/,  or 
lliat  the  Sj  are  independent  of  the  Examples  which  illustrate  the  applications  of  the  theory 
are  included. 

(3)  Consider  a  set  of  components  which  usually  work  together  as  pans  in  some  complex  system. 
It  is  often  of  interest  to  study  the  behaviour  of  one  of  the  components  (or  of  a  subset  of  the 
components)  “separately”  from  tlie  other  components.  This  leads  one  to  consider  die  marginal 
distribution  of  the  panicular  component  or  of  the  set  of  components.  However,  formally,  the 
classical  definition  of  marginal  distributions  is  not  always  meaningful  in  some  applications  in 
reliability  dieory.  This  difficulty  with  the  classical  definition  of  marginal  distribution  led  us  to 
study  a  new  notion  of  marginal  distributions.  The  new  notion  is  based  on  die  notion  of  multi¬ 
variate  imperfect  repair  which  the  Pi’s  have  developed  earlier.  The  study  is  described  in  [7], 
The  study  shows  how  to  formalize  different  interpretations  of  “separability”  of  components. 
In  particular,  the  study  also  shows  how  to  model  conjunctions  of  reliability  systems  (that  is, 
how  to  have  a  formal  way  of  modeling  combinations  of  components  from  different  systems  at 
some  point  in  time  in  which  tliey  are  not  new  anymore  and  may  have  different  ages). 

(4)  The  use  of  the  notion  of  multivariate  minimal  repahs  has  already  been  mentioned  above  as  a 
tool  for  modeling  dynamic  conditional  marginal  distributions.  Minimal  repairs  are  also  of  in¬ 
terest  in  situations  of  tight  budgets  which  do  not  permit  a  perfect  replacement  each  lime  a  com¬ 
ponent  fails.  In  such  case  the  component  is  repaii’cd  so  that  it  is  brought  back  to  a  state  of 
being  functional,  but,  in  general,  it  is  not  as  good  as  new  then.  In  some  such  situations  the 
budget  may  not  even  permit  an  unlimited  number  of  minimal  repairs.  It  may  suffice  only  for 
some  fixed  number  of  minimal  repairs  which  in  such  a  case  need  to  be  allocated  for  the  vari¬ 
ous  components  of  the  system.  The  question  of  what  i.s  the  optimal  allocation  then  arises.  In 
[12]  we  have  studied  questions  of  this  kind,  and  we  have  found  the  optimal  allocations  in  a 
variety  of  tight  budget  situations.  Not  only  minimal  repairs  have  been  considered,  but  also  the 
optimal  policies  of  assigning  various  kinds  of  standby  replacements.  The  notions  of  minimal 
and  imperfect  repairs  arc  u.seful  as  well  in  other  recent  studies  done  under  tins  tliis  Grant.  For 
example,  in  [5]  these  notions  arc  used  as  a  model  for  repair  of  discrete  time  components.  In 
fact,  one  can  view  the  results  in  [5)  as  a  po.s.siblc  version  of  imperfect  repair  in  di.screte  time. 

Some  work  that  has  not  been  described  in  the  proposal  has  also  been  done  under  this 
Grant.  In  [13],  [14]  and  [15]  various  notions  of  stochastic  convexity  and  stochastic  majorization 
are  studied  in  detail.  In  [16]  and  [17]  it  is  shown  that  often  (but  not  always)  tlie  statistical  informa¬ 
tion  content  of  positively  dependent  random  variables  is  larger  ilian  the  statistical  information  con¬ 
tent  of  comparable  independent  random  variables. 


0CM9-93  TUE  10:52 


P.05 


REFERENCES 

[1]  Shaked,  M.  and  Shanthikumar,  J.  G.  (1991),  Dynamic  multivariate  mean  residual  life  func¬ 
tions,  Journal  of  Applied  Probability,  28,  613-629. 

[2J  Shaked,  M.  and  Shanthikumar,  J.  G.  (1993),  Multivariate  conditional  hazard  rate  and  mean 
residue  life  functions  and  thtii  ippWciUonfi,  Reliability  and  Decision  Making  (Eds:  R.  E. 
Barlow,  C.  A.  Clarotti  and  F.  Spizzichino),  137-155. 

[3]  Shaked,  M.  and  Shanthikumar,  J.  G.  (1991),  Dynamic  construction  and  simulation  of  random 
vectors,  Topics  in  Statistical  Dependence  (Ed.s:  H.  W.  Block,  A.  R.  Sampson  and  T.  H. 
Savits),  415-433. 

[4]  Shaked,  M.,  Shanthikumar,  J.  G.  and  Valdcz-Tonres,  J.  B.  (1993),  Discrete  hazard  rate  func¬ 
tions,  Computers  and  Operations  Research,  to  appear. 

[5]  Rocha-Martinez,  J,  M.  and  Shaked,  M.  (1992),  A  discrete- time  model  of  failure  and  repair, 
Technical  Report,  Deparunent  of  Mathematics,  University  of  Arizona. 

[6]  Shaked,  M.,  Shanthikumar,  J.  G.  and  Zhu,  H.  (1993),  Discrete  repairable  reliability  systems. 
Naval  Research  Logistics,  40,  769-780. 

[7]  Shaked,  M.  and  Shantliikumar,  J.  G.  (1993),  Dynamic  conditional  marginal  distributions  in 
reliability  theoiy.  Journal  of  Applied  Probability,  30,  421-428. 

[8]  Shaked,  M.,  Shantliikumar.  J.  G.  and  Valdcz-Torres,  J.  B.  (1991),  Discrete  probabilistic  or¬ 
derings  in  reliability  theory.  Technical  Report,  Department  of  Mathematics,  University  of 
Arizona. 

[9]  Shaked,  M.  and  Zhu,  H.  (1992),  Some  results  on  block  replacement  policies  and  renewal  the¬ 
ory,  Journal  of  Applied  Probability,  2  9,  932-946. 

[lOJPellercy,  F.  and  Shaked,  M.  (1993),  Stochastic  comparison  of  some  wear  proces.scs. 
Probability  in  the  Engineering  and  Informational  Sciences,  7,  421-435. 

[llj  Shaked,  M.  and  Shanthikumar,  J.  G.  (1991),  Dynamic  multivariate  aging  notions  in  reliability 
theory.  Stochastic  Processes  atui  Their  Applications,  38,  85-97. 


-4- 


0CM9-93  TUE  10:52 


P.  06 


[12J  Shaked,  M.  and  Shanthikumar,  J.  G.  (1992).  Optimal  allocation  of  resources  to  nodes  of  par¬ 
allel  and  series  systems,  Advo/icer  in  Applied  Probability ,  24,  894-914. 

[13]  Li,  H.  and  Shaked,  M.  (1993),  Stochastic  majorization  of  stochastically  monotone  families  of 
random  variables.  Advances  of  Applied  Probability,  to  appear, 

[14]  Li,  H.  and  Shaked,  M.  (1993),  Stochastic  convexity  and  concavity  of  Markov  processes. 
Mathematics  of  Operations  Research,  to  appear. 

[15]  Shaked,  M.  and  Shanthikumar,  J.  G.  (1991),  Regular,  sample  path  and  strong  stochastic  con¬ 
vexity:  A  review.  Stochastic  Orders  and  Decision  under  Risk  (Eds:  K.  Mdsler  and  M. 
Scarsini),  320-333. 

[16]  Shaked,  M.  and  Tong,  Y.  L.  (1992),  Comparison  of  experiments  via  dependence  of  normal 
variables  with  a  common  marginal  disuibution.  Annals  of  Statistics,  20,  614-618. 

[17]  Shaked,  M.  and  Tong,  Y  L.  (1992),  Comparison  of  experiments  of  some  multivariate  distri¬ 
butions  with  a  common  marginal.  Stochastic  Inequalities  (Eds:  M.  Shaked  and  Y.  L.  Tong), 
388-398. 


-5- 


